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ABSTRACT 



r \ ' The outer haloes of the Milky Way (MW) and Andromeda (M31) galaxies contain 

. I as much important information on their assembly and formation history as the prop- 

i-Q ■ erties of the discs resident in their centres. Whereas the structure of dark matter (DM) 

Mh' haloes has been studied for a long time, new observations of faint structures hiding 

Q , in the depths of the stellar halo have opened up the question of how the stellar halo 

$— ( ■ is related to the DM underlying it. In this paper we have used the Constrained Local 

(/3 I UniversE Simulation (CLUES) project to disentangle the stellar and DM component 

C^ . of three galaxies that resemble the MW, M31 and M33 using both DM only simula- 

tions and DM + gas-dynamical ones. We find that stars accreted in substructures and 
then stripped follow a completely different radial distribution than the stripped DM: 
^ I the stellar halo is much more centrally concentrated than DM. In order to understand 

\jQ ' how the same physical process - tidal stripping - can lead to different ^ = radial pro- 

V^ i files, we examined the potential at accretion of each stripped particle. We found that 

CO ■ star particles sit at systematically higher potentials than DM, making them harder 

■^ ' to strip. We then searched for a threshold in the potential of accreted particles 0th, 

r — , above which DM particles in a DM only simulation behave as star particles in the 

f^ ' gas-dynamical one. We found that in order to reproduce the radial distribution of star 

particles, one must choose DM particles whose potential at accretion is > 16(/)subhaio, 
where ^subhaio is the potential at a subhaloes edge at the time of accretion. A rule 
as simple as selecting particles according to their potential at accretion is able to re- 
produce the effect that the complicated physics of star formation has on the stellar 
distribution. This result is universal for the three haloes studied here and reproduces 
^ ' the stellar halo to an accuracy of within ^2%. Studies which make use of DM parti- 

cles as a proxy for stars will undoubtedly miscalculate their proper radial distribution 
and structure unless particles are selected according to their potential at accretion. 
Furthermore, we have examined the time it takes to strip a given star or DM particle 
after accretion. We find that, owing to their higher binding energies, stars take longer 
to be stripped than DM. The stripped DM halo is thus considerably older than the 
stripped stellar halo. 

1 INTRODUCTION bottom up l| White fc Ree^l 19781 : [Pwis et al.|[l985l ) - the first 

objects to collapse at high redshift are small sub-galactic 

, ,. , ,, , , , , . , , , ,, units, that condense out of small perturbations in the initial 

According to the current accepted cosmological model, the , . r. n n mi i i • i • i • i 

, p „„rw , , , , ,, /T^»r\ r^niw density held, these clumps then merge m an hierarchical 

universe is composed ot 26% cold dark matter (DM), 70% r -, ■ , -, , i , . , 

, , ,,, , .(w , / [t; ', — , — nUr.^rJx tashion to construct the large bound obiects we observe m 

dark energy (A) and 4% baryons (e.g. ISpergel et al.M2007l ). , , , . ,,..,,,. 

„, , . ,\ ,, , » ^T^, , , r r ,1 the local universe. Visible luminous matter - stars - are be- 

atructure m the so-called ACDM cosmology forms from the 

© 0000 RAS 



X 



C^ 



Libeskind et al. 



lieved to be formed when Giant Molecular Clouds collapse 
in the potential wells of these bound blobs of DM and gas. 

The merging process gives rise to DM haloes, which 
today host bright central galaxies such as the Milky Way 
(MW) and the Andromeda galaxy (M31) in their cores. The 
outskirts of such DM haloes are populated by a two compo- 
nent medium: diffuse matter and matter bound to substruc- 
tures. Much of the mass is found bound to satellite galaxies 
which orbit within their parent halo. The properties (age, 
orbital parameters, spatial distribution, kinematics, etc) of 
luminous satellite galaxies can teach us a lot regarding the 
formation of their hosts and have been the target of nu- 
merous observational and theoretical studies. Indeed, the 
past 5 years has seen an increased focus on the detection 
of satellite galaxies and has resulted in around a dozen new 
satellite s being detect ed by the SDSS JKoposov et al.ll20o4 
[Belo ku rov et all l2009t ) . 

Yet the z = Q satellite galaxy population is not a full 
survey of all the substructures accreted by the parent DM 
halo, since many substructures accreted at high redshift will, 
by 2: = 0, have been tidally disrupted by the host poten- 
tial, resulting in the stripping of dark matters and stars. 
Indeed it is believed that the stellar halo - stars exterior 
to the central galaxy and not bound to substructures - was 
for med by the tea r ing of stars from accreted satellite galax- 
ies. ICooper et al.l ((20101) argue that the vast majority of 
stars in the M W's halo were stripp ed from just one or two 
large satellites. IZolotov et al.l (|2009|) have studied the stellar 
halo in gas-dynamical/N-body simulations and have iden- 
tified that in fact the stellar halo has a dual-origin: part 
of it was created via tidal stripping of stars from disrupted 
satellites, and part was pushed out of c entral galaxies during 
minor mergers. In a follow up paper, IZolotov et al.l (|2010l ) 
argue that the metallic abundance patterns (of [Fe/H] and 
[0/Fe]) of stars can be used to distinguish between theses 
different formation mechanisms. 

Observations using the SDSS by, e.g.. lBell et"aLl (|2008l ') 
have indicated that the MW's halo is consistent with be- 
ing formed entirely out of accreted de bris material. This is 
in disagreement with observations by ICaroUo et al.l (|2008l ) 
who find clear differences between kinematical properties of 
the inner and outer stellar halo - stars in the inner halo are 
found to exhibit a net prograde rotation while the outer halo 
is dominated by retrograde motion. It has been suggested 
that the differences in net rotation of stars in the halo be- 
tray a dual origin, a result recently supported bv lBeers et al.l 
l|201ll ) who find differences in kinematics for inner and outer 
halo stars. 

The situation with M31 is similar as observations seem 
to favor an accreted or i gin. B y studying the ages of stars in 
the halo, iBrown et al.l (|2008l ) argue against a in-situ origin 
and, because the stars are by and large relatively old - they 
argue for an hiera rchical build up of M 31. By focusing on the 
metal enrichment. [Gilbert et al.l (|2009l ') too seem to argue for 



an entirely accreted stellar halo with little evidence of in-situ 
star formation. 

The assembly history of the DM halo on the other hand, 
is more difflcult to pin down as direct observ ations are by 



defin i tion impossible. Yet many authors (e . g . Klvpin et al 



200 ll: IWechsler et~aIll20Q2l : IZhao et al.ll2003l : IPiemand etal 
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20071 . and references therein) have used A'^-body simula- 
tions to determine the relative importance of the two main 
"modes" of h alo growth: d i ffuse accretion versus mergers. 
Most recently IWang et al.l (120101 ) have used the Aquarius 
simulation and found that ambient accretion contributes the 
largest amount of material to the dark halo. 

It is thus unclear if the conclusions drawn from obser- 
vations of our stellar halo - that the stars were stripped 
from infalling satellite galaxies - are consistent with DM 
only simulations which point towards a diffuse smoothly ac- 
creted halo where mergers and debris material play a minor 
role. In this paper we disentangle these two components in 
order to understand how they co-evolved. 



2 METHODS 

In this section we describe in brief the simulations used as 
well as the halo and subhalo finding algorithm employed to 
identify satellites. 



2.1 Constrained Simulations of the Local Group 

The simulations used in this work are embedded in 
the Constrained Local UniversE Simulation (CLUES) 
project and have be en already studied i n a number 
of re c ent papers (e.g Libeskind et al. 2010l: Knebe et al.l 



2OIOI; iKlimentowski et al.l I2OI0I: iLibeskind et al.1 120111 



Knebe et al.ll201l | ) We refer the reader to those papers (in 



particular ILibeskind et al.ll2010l ) for details on how the con- 
straints were generated and how the simulations were run: 
we highlight just the salient points here for clarity. 

We choose to run our simulations using standard ACDM 
initial conditions, that assume a Wilkinson Micr owave 
Anisotropy Probe 3 cosmology (|Spergel et al.l |2007| ) , i.e. 
Q.m = 0.24, fib = 0.042, ^A = 0.76 and h = 0.73. We use 
a normalization of erg = 0.75 and an n = 0.95 slo pe of the 
powe r spectrum. We use the MPI code GADGET2 (jSpringell 
[2OO5I) to simulate the evolution of a cosmological box with 
side length of Lbox = 64/i,~ ^Mpc, befor e applying the zoom 
technique (see e.g. iKlvpin et al.ll200ll l around a region of 
interest. 

Instead of seeding our initial conditions as a just a ran- 
dom cube of space, the initial conditions of our volume are 
constrained to reproduce, at z = a number of objects that 
compose the local environment (see lHoffman fc Ribaklll99ll . 
for details on how constraints initial conditions are gener- 
ated), including a "virgo" cluster, a "coma" cluster and a 
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"Local Group" . Our method allows us to properly constrain 
the large scales (i.e. those still linear by 2: = 0) but we do 
not constrain the local group itself. In order to obtain a lo- 
cal group in the correct environment, three low resolutions 
constrained simulations are run with varying random seeds. 
Each z = low resolution simulations is then examined 
and if an object that resembles the local group is found (by 
construction this will be in the correct place), these initial 
conditions are selected for high resolution re-simulation. 

Our initial density field includes both DM and gas par- 
ticles. Under the right conditions, gas particles may spawn 
star particles which interact gravitationally in the same 
way as DM (i.e. as point particles with a given softening 
length). Each gas particle may have up to two star forma- 
tion episodes, each time spawning a star of half its original 
mass. In order to conserve mass we reduce the gas parti- 
cle's mass each time a star particle is spawned, resulting in 
gas particles that have one or two times the mass of star 
particles (corresponding to gas particles that have spawned 
one or no star particles). When a gas particle spawns its 
second star particle, it ceases to exits. Star particles repre- 
sent stellar populations and are given the metalicity of the 
gas particle that spawned it. The massive stars born (with 
M > IOMq) in this population explode instantaneously as 
supernovae type II, polluting the environment with metals 
and producing stellar winds. More details on the star for - 
mation prescription can be found in (jLibeskind et al.ll2010r ). 

We resimulate just the region of interest around the lo- 
cal group. We centre a sphere of radius 2 h~^Mpc around the 
local group and populate it with ~ 5.2 x lO'' low mass, high 
resolution particles. Within our local group we are thus able 
to achieve a particle mass of just Afdm ~ 2.54 x 10^ h~^MQ 
for DM and Af^tar = 2.21 x lO"* /i"^Mq for star particles. 

Our constraints reproduce a cosmography which closely 
resembles the observed Local Group. In Table. [T]we compare 
properties of the simulated local group with observations of 
the real onqj. Although our results do not match the ob- 
servations perfectly, the cosmography simulated using our 
constraints captures the essence - in terms of mass and dis- 
tances - of the observed Local Group. 

In addition to our gas dynamical SPH simulation, we 
also have a DM only version seeded from the same ini- 
tial conditions. A comparison b etween the two s imula tions 
has already b een highlighted in lLibeskind et al.l (|2010|) and 
iKnebe et al.l ((2010). The DM only simulation has similar 
spatial and mass resolution and reproduces the same three 
main haloes as the gas-dynamical simulation. We use the 
DM only simulation solely in Section 13.21 where we try to 
find a recipe by which particles in a DM only simulation can 
be used to reproduce the radial distribution of the stellar 



^ In a future paper we intend to study in detail the cosmography 
produced by our constrained simulations. 



Property 


Simulated LG 


Observed LG 


Reference 


A^MW 


6.57 X 1O"M0 


IO12 Mq 


[1,2,3] 


Mm31 


8.17 X 10" M0 


8.2 X IO^Mq 


[4] 


Mm33 


2.02 X 10^1 Mq 


6 X 10^0 AfQ 


[5] 


ryrw 


220 kpc 


253 kpc 




^M31 


245 kpc 


237 kpc 




^M33 


183 kpc 


100 kpc 





Table 1. The z = properties of the simulated and observed 
Local Group. From the top row own, we show the following prop- 
erties: the mass of the MW's halo (Afiviw), the mass of M31's halo 
(Aflviai), the mass of M33's halo (A/M33), the virial radius of the 
MW halo (rMw)i the virial radius of M31's halo (rMsi), and the 
virial radius of M33's halo (^Mss)- Note that the "observed" virial 
radii are calculated from the observed virial m asses; they are thu s 
unref erenced. The refer ence s are as follows: [1] Xue et al.l ll2008r: 
[2] Klvpiri et al.l ||2002|'I: [3] ISmith et all ||2007| '): [4] ISeigar et al.] 
Il2008l) :f5l ICorbeUil 1I200I ') 



halo, without the necessity of a complicated semi-analytical 
model. 



2.2 The halo and subhalo finding algorithm 

In this section, we explain how our halo and subhalo find- 
ing algorithm works. In order to identify haloes and sub- 
haloes in our simulation we have run the MPI-f OpenMP 
hybrid halo finder AHF (AMIGA halo finder, to be downloaded 
freely fro m http : //popia . ft . u am . es /AMIGA) described in 
detail in iKnoUmann fc Knebj (120091). AHF is an improve- 
ment of the MHF halo finder (|Gill et aLll2004l L which locates 
local over-densities in an adaptively smoothed density field 
as prospective halo centres. The local potential minima are 
computed for each of these density peaks and the gravita- 
tionally bound particles are determined. Only peaks with at 
least 20 bound particles are considered as haloes and re- 
tained for further analysis. In practice for this work, we 
only consider subhaloes with more than 100 particles. We 
would like to stress that our halo finding algorithm auto- 
matically identifies haloes, sub-haloes, sub-subhaloes, etc. 
For more details on the mode of operation and actual func- 
tionality we refer the reade r to the code description paper 
(|Knollmann fc KnebdlJOOgl L 

For each halo, we compute the virial radius rvir, that is 
the radius r at which the density M{< r)/(47rr"^/3) drops 
below AvirPback. Here pback is the cosmological background 
matter density. The threshold Avir is computed using the 
spherical top-hat collapse model and is a function of both 
cosmological model and time. For the cosmology that we are 
using, Avir = 355 a.t z = 0. 

Subhaloes are defined as haloes which lie within the 
virial radius of a more massive halo, the so-called host halo. 
As subhaloes are embedded within the density of their re- 
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spective host halo, their own density profile usually shows a 
characteristic upturn at a radius rt < rvir, where rvir would 
be their actual (virial) radius if they were found in isolationU 
We use this "truncation radius" rt as the outer edge of the 
subhalo and hence subhalo properties (i.e. mass, density pro- 
file, velocity dispersion, rotation curve) are calculated using 
the gravitationally bound particles inside the truncation ra- 
dius rt- For a host halo we calculate properties using the 
virial radius r-vir- 

We build merger trees by cross-correlating haloes in 
consecutive simulation outputs. For this purpose, we use 
a tool that comes with the AHF package and is called 
MergerTree. As the name suggests, it serves the purpose 
of identifying corresponding objects in the same simulation 
at different redshifts. We follow each halo (either host or 
subhalo) identified at redshift z = backwards in time, 
identifying as the main progenitor (at the previous redshift) 
the halo that both shares the most particles with the present 
halo and is closest in mass. The latter criterion is important 
for subhaloes given that all their particles are also typically 
bound to the host halo, which is typically orders of magni- 
tude more massive. Given the capabilities of our halo finder 
AHF and the appropriate construction of a merger tree, sub- 
haloes will be followed correctly along their orbits within the 
environment of their respective host until the point where 
they either are tidally destroyed or directly merge with the 
host. 



"in-situ"; if it formed in any other halo we say it formed 
"ex-situ" . 

Ex-situ stellar halo particles have thus been stripped 
from the accreted substructure within which they were born 
and orbit within the parent as debris material. We thus also 
refer to these particles as "stripped" star particles and use 
this term interchangeably with "ex-situ". In situ halo star 
particles were formed in the centre of the progenitor of the 
host halo and are then pushed out by merging or migrational 
processes fe.g. IZolotov et al.ll2009l '). 

For the diffuse DM halo, particles may also be stripped 
from accreted subhaloes. Unlike star particles however DM 
particles may also be "smoothly" accreted, in other words 
accreted by the main progenitor either individually from the 
ambient cosmic background or in substructures below our 
subhalo resolution limit of 100 particles. 

In this work we have used the term stripped particles 
to denote those particles that become unbound from the 
subhalo in which they were accreted by any physical mech- 
anism. Our term is general in that we do not differentiat e 
between tidal or resonant stripping (jD'Onghia et al.ll2009l ). 

Throughout this paper we choose not to stack our three 
haloes into single plots, since the variety of results obtained 
is significant and due to the different and unique merger 
histories of each halo. 



2.3 Identifying the stellar halo 

In this section we explain the nomenclature used for the 
different particle sets in our analysis. We exclude from our 
analysis the inner baryonic component (i.e. the galactic disc) 
and use the term "outer halo" to refer to the region between 
0.1 rvir and rvir. We identifying those DM and stellar parti- 
cles that at z = are within this region and then excise all 
particles bound to substructures. We call the "swiss cheese" 
like remains the diffuse (stellar or DM) halo. We focus our 
analysis on the origin of this diffuse component. 

We refer to star particles in the diffuse outer halo as the 
"stellar halo" . For each of these particles, our simulation pro- 
vides us with the age of the universe when it formed. Since 
star particles by construction can only be formed in high 
density environments (and are thus bound to a subhalo at 
the moment of their birth) we can use the star particle's age 
to locate the appropriate snapshot and thus the (sub)halo in 
which the star was spawned. If the star formed in the main 
progenitor of its 2 = host, we say the particle formed 

^ Please note that the actual density profile of subhaloes after the 
removal of the host's background d rops faster than for isolated 
haloes (e.g. lKazantzidis et al.ll2004r) : only when measured within 
the background still present will we find the characteristic upturn 
used here to define the truncation radius rt- 



3 RESULTS 

We begin by studying the radial distribution of mass within 
each (DM and stellar) component. In Fig. [T] we show the 
diffuse mass interior to a given radius as a function of radial 
distance from the centre, for DM (black solid) and stars 
(red solid) for our three galaxy haloes. We normalise each 
component by the total diffuse halo mass within the outer 
halo. Note that the diffuse DM makes up roughly 80% of 
the halo's full virial mass (jGao et al.ll2004l) and i s used in 
the semi-analytical investigations of ICooper et al.l ((20101) as 
a proxy with which to study the stellar halo. 

When comparing the stellar component to the DM com- 
ponent (the solid red line to the solid black line), a very 
stark difference is immediately visible. Although star and 
DM particles are treated equally in the simulations gravity 
calculation, their z = distribution differs dramatically in 
that the star particles are highly concentrated towards the 
centre of the outer halo while the DM roughl y follows an 
NFW distribution (not shown here, but see e.g. IWang et al.l 
l2010t ). 

The extreme central concentration of stars is likely due 
to radiative cooling, which causes gas particles to lose en- 
ergy, fall to the centre of the halo where the densities are 
high enough for star formation, adiabatically contract the 
DM and dee pen the potential, a resu lt that has been know 
since at least iBlumenthal et al.l (|l986l ). Yet star particles can 
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Figure 1. The cumulative mass profile (i.e. total mass fraction within a radius r) for the diffuse component of our MW (upper left), 
M31 (lower left) and M33 (upper right) outer haloes. Both DM (black) and stellar (red) curves are normalized by the total mass in the 
respective diffuse component within the outer halo. The solid red curve shows the mass profile of the diffuse stellar component. The red 
dot-dashed line shows the mass profile for those stars born within the main progenitor, known as "in-situ" stars. The dashed red line 
shows the mass profile for stars accreted in clumps and later stripped from them such that at 2 = they are bound just to the main halo. 
Similarly, the solid black curve shows the mass profile of the diffuse DM component. The dot dashed black curve shows this quantity for 
smoothly accreted DM and the dashed black line represents the mass profile for stripped debris material. 
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also fall to the centre by losing angular momentum through 
dynamical friction against the halo background (as can the 
DM) a process kn own to be more effe c tive for clumps that 
contain baryons (JEl-Zant et alj I2OOII : iRomano-Diaz et alj 
l200a ). These two routes of galaxy formation are both impor- 
tant a nd have been studied in great detail by e.g. lOser et al.l 
(|2010l ). Thus in order to better understand how this differ- 
ence in radial profile of the diffuse component has arisen, we 
must examine in greater detail each component individually. 

We now focus on the stellar component, specifically on 
its two constituent subsets: in-situ and ex-situ star particles. 
The stripped and in-situ stars display considerably different 
radial profiles (red dotted and red dashed curves respec- 
tively in Fig. [TJ , perhaps not surprising given their different 
physical origins. Whereas the cumulative in-situ stellar mass 
rarely grows beyond 0.2 rvir, the stripped stars are found in 
abundance at all radii throughout the stellar halo. 

Note that in-situ stars contribute very little to the total 
mass when averaged over the entire diffuse halo. The outer 
haloes are dominated by ex-situ star particles that comprise 
- 80% of the MW's outer halo and ^ 90% of M3f 's and 
M33's outer halo. The exact amount is due to the unique 
merger history of each individual halo. Because the outer 
halo is mostly made up of these stripped particles, their 
mass profiles dominate the total mass profile of the outer 
halo. The fact that the stellar halo is composed primarily 
of stripped stars bodes well since this population has direct 
counterpart in the DM component. 

We now look at the two components of the diffuse DM 
halo, the smoothly accreted and stripped particles. Note 
that the smoothly accreted component has no stellar coun- 
terpart in since a negligible number of star particles are "pre- 
stripped" and smoothly accreted. From Fig. [T] we see that 
the smoothly accreted DM constitutes the major part of the 
outer diffuse DM halo, contributing ~ 60 — 70% to its mass. 
Furthermore, its radial distribution has roughly the same 
shape as the total outer DM halo (not surprisingly since it 
dominates the halo's mass) and is markedly different from 
the stripped DM debris. 

Unlike the smoothly accreted DM, the stripped DM 
(shown in Fig. [T] as the dashed black line), has a direct 
counterpart among the star particles (see above). Yet the 
two stripped profiles have completely divergent shapes and 
it is difficult to find any similarities between the two popu- 
lations. Furthermore, stripped DM contributes just ~ 40%, 
~ 30% and ~ 25% to the total mass of the outer halo for 
the MW, M3f and M33 respectively (versus ~ 80 - 90% for 
stripped stars). 

But how different are the radial profiles of just the 
stripped stellar and DM components when normalized to 
their respective rvir masses? In Fig. [2] we plot the radial 
profiles of these two components. The two curves clearly de- 
viate substantially from each other, both qualitatively and 
quantitatively. Stripped stars (red lines of Fig. [2| are, as 



noted earlier, significantly more centrally concentrated then 
the stripped DM (black line). Note that the cause for this 
discrepancy cannot be the deeper potential of the host - adi- 
abatic contraction of the DM due to radiative cooling of the 
gas only affects the region interior to O.Irvir (see Fig. 2 of 
iLibeskind et aLlbOloD 

Thus despite the fact that these two populations are 
born out of the same physical processes (and treated equally 
by the simulations gravity calculator), they still have at 
z = G very dierent radial distributions. We now examine 
the origin for this dichotomy. 



3.1 Examination of the potential at accretion 

The chance that a particle will at some point become 
stripped from the clump that it was bound to at accretion, 
depends on a variety of properties of the sub-clump (for ex- 
ample its orbit, internal structure, spin, etc). One important 
factor is the nature of the potential well of the subhalo that 
each accreted particle sits in at accretion. The potential of 
each particle at accretion can be calculated by assuming 
that the accreted subhalo obeys spherical symmetry. The 
potential at a distance r from the subhalo's centre is thus: 



= G 



M{< 



-dr' + 0(0) 



(1) 



where 0(0) is normalized s uch that the potential is nu ll at 
infinity (cf the appendix of iKnoUmann fc Knebell2009l ). At 
the subhaloes edge the potential is 



GMu 



(2) 



Since the depth of a potential well at the centre of a 
subhalo is a measure of the concentration (which in turn 
depends on other subhalo dependent pro perties like the to- 
tal mass or formation tim e - see e.g. IZhao et al.l l2003l : 
iMunoz-Cuartas et al.ll201ll . among others) we normalize the 
potential of each particle by the potential at the subhaloes 
edge to obtain a dimensionless number that is more or less 
independent of global subhalo properties. 

We plot the distribution of the logarithm of this value 
in Fig. [3]for stars (red) and DM (black). The vertical dashed 
lines indicate the median values, which for the stars are 
always more than fOX(j!>aubhaio- This means that the stars 
that get stripped and which subsequently end up dominat- 
ing the stellar halo had, at the time of accretion, a potential 
ten times greater then that at the edge of their host sub- 
halo. This is a reflection of the fact that stars are "more 
bound" than other particles in the same subhalo; yet there 
is some tension because the more bound a particle is to 
its subhalo, the less likely it is to be stripped from it. The 
stripped DM, on the other hand, occupies "less bound" re- 
gions of the subhalo's potential at accretion, evident by the 
fact that the median potential is always significantly lower 
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Figure 2. The cumulative mass profile (i.e. total mass fraction within a radius r) for just the stripped component of the diffuse haloes of 
our MW (left), M31 (centre) and M33 (right) haloes. Both DM (black) and stellar (red) curves are normalized by the total mass in the 
respective stripped component within r^ir. The dashed blue curve corresponds to those DM particles in the gas-dynamical simulation 
whose potential at accretion (t^sph) was greater than 10^'^^®, lO^*^^®, 10^"^'-* times that of the host (^subhalo) they were accreted in. 
The dashed black line corresponds to those DM particles in the DM only simulation whose potential at accretion (</>dm) was greater 
than lO^-^^^, lO^'l^^, and IQl-^SS times that of the host for the MW, M31, and M33 respectively. See section |32] for more on how these 
thresholds were obtained. 



than that for the star particles. It may be counter-intuitive 
that so many stripped particles have apparently high values 
of log (/)acc/0aubhaio, but this is a direct result of the fact that 
particles in subhaloes are not uniformly distributed; instead 
the majority of the particles that make up a subhalo are con- 
centrated in the centre of the subhalo, where the potential 
is much higher than at the subhalo's edge. 

In essence, Fig. [3] reflects the distribution of binding 
energies at accretion of stripped particles. We could have 
plotted this value instead, but by normalizing by the poten- 
tial at the halo's edge we obtain a dimensionless quantity 
which is subhalo independent. 

Fig. [3] reveals that DM particles are stripped from the 
fluffier outer regions of a subhalo, while star particles are 
stripped from the deep interior of their subhaloes. This can 
explain why the number of stripped star and DM particles is 
so different: not only are there fewer stars to strip to begin 
with, but they are harder to strip after accretion. 



3.2 Finding a DM subset that mimics the stellar 
halo 

We wish to thus identify a subset of the DM that follows at 
z — the same radial distribution as the stripped stars but 
that does not - in any way - depend on stellar properties. 
In this way, pure DM-only simulations may be used and 
a subsample that accurately reproduces the stellar radial 
distribution can be obtained. 

From Fig.[3]we know that stars were more tightly bound 
in bigger subhaloes at the moment of accretion, while the 
stripped DM was more loosely bound to smaller subhaloes 



when they fell in. We thus wish to choose DM particles 
that sit at the same depth (or deeper) of the potential 
well as the stripped stars. In principle we could easily se- 
lect a sub-sample of DM particles whose distribution of 
log(^acc/'^subhaio perfectly matches that of the stars; yet we 
wish our criteria for sub-selecting the DM to be independent 
of any stellar properties. 

We thus define a threshold (j>th for the value of 
<^acc/<^subhaio above which all DM particles are selected as a 
proxy for stars. After defining a threshold and obtaining a 
DM sub-sample, we examine the radial distribution of these 
particles and compare it (by calculating x^) to the star par- 
ticle curves in Fig. [21 By smoothly varying (f>th and exam- 
ining the value of x^ for each halo, we are able to obtain a 
DM sub-sample which very closely matches the stellar radial 
distribution. 

Although the optimum (fith varies slightly across our 
three haloes, it does so weakly having values of i^th mw = 
lO^-^®^ « 15, <t>tb M31 = 10'°^® « 11, and (^th m33 = 
10^°™ « 12, for the MW, M31 and M33 respectively. We 
plot the radial distribution of DM particles that meet the 
(/)th criteria in Fig. [2] as the dotted blue line, and note that 
by selecting these subsets we nearly perfectly recover the 
radial stellar distribution for each halo. Note that roughly 
the same number of DM particles meet these criteria as halo 
star particles. 

Yet in order to develop a rule by which DM only simula- 
tions can be used to study the stellar halo, we must attempt 
to match the radial distribution of a subset of DM particles 
in a DM only simulation to the stripped star particles in the 
gas-dynamical one. The thresholds mentioned above work 
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Figure 3. The cumulative (top) and differential (bottom) distribution of the logarithm of a particle's potential at accretion (ipacc) 
normalized to the value of the potential at the subhalo's virial radius ((/'subhalo)- In black we show this distribution for DM and in red 
for star particles. The median values are shown by the vertical dashed lines and arc lO^'""^ (lO"'®^), 10^'^^ (lO"'*^), and lO"'^'^^ (lO"'^'^), 
for the stars (DM) in the MW, M31 and M33 haloes. 



well for the gas-dynamical simulation but, since the poten- 
tial of each halo is aected by adiabatic contraction due to 
the collapsed baryons, the limit is not directly comparable 
to DM only simulations. 

In order to address this concern, we perform the exact 
same analysis we have thus far presented on our DM only 
simulation. We perform the same x^ minimization test to 
obtain the best fit radial profile of DM particles in the DM 
only simulation. We show this distribution as the dashed 
black line in Fig. (2] In order to reproduce the stripped stellar 
distribution we need to select DM particles that are slightly 
deeper than their counterparts in the gas dynamical simula- 
tion, having potentials greater than cj>th mw = 10^'^^^ « 17, 
0th M31 = 10'-''' ^ 14, and <?!.th m33 = lO'''"' ^ 20, for the 
MW, M31 and M33 respectively. The necessity of select- 
ing particles at a deeper part of the potential in DM only 
simulations is due to the relatively shallower potentials in 
simulations without baryons. 

We note that the number of stripped DM halo particles 



that meet the (jisubhaio criteria in the gas dynamical simula- 
tion is: ~ 2%, 5% and 6% of the total stripped halo for the 
MW, M31 and M33 respectively. This fraction is roughly the 
same as is found in the DM only simulations where ~ 1%, 
3% and 3% of the diffuse stripped halo meet the criteria for 
the MW, M31, and M33 respectively. Furthermore the abso- 
lute number of particles that meet these criteria is roughly 
the same across all three haloes and is the same order of 
magnitude as the number of stars in the stripped stellar 
halo. 

We now attempt to obtain a universal threshold in the 
potential at accretion of DM particles in DM only simula- 
tions that reproduce the stripped stellar profile at z = 0. We 
begin by averaging the stripped stellar profile of our three 
haloes and plotting it as the red line in Fig. U) We then co- 
add the diffuse stripped DM halo particles and select DM 
particles according to whether they are above or below the 
threshold potential (/)th at accretion. We examine three fidu- 
cial values for the threshold: 0th = 150subhaio, 160subhaio 
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Figure 4. The co-added eumulative diffuse halo radial profile at 
2 = averaged over our 3 haloes. In red we show the stellar distri- 
bution. The black lines indicate different thresholds in the poten- 
tial, above which DM particles in our DM only simulation were 
selected. The dashed line is thus the radila distribution of DM 
particles who at accretion had potentials greater than 15</>subhalo 
times the potential at the subhaloes edge. The dot - dashed line 
is for potentials greater than 16</<subhalo and the triple dot dashed 
line is for potentials greater than 17</<Bubhalo- The bottom panels 
shows the residuals when the stellar distribution is subtracted. 



and 170subhaio, and plot the radial distribution of DM par- 
ticles that meet these criteria as the dashed, dot dashed and 
triple dot dashed lines in Fig. [l] In the bottom panel of 
Fig. 3] we show the residual value when the DM profile is 
subtracted from the stellar profile. We note that the best 
value is (f>th = 16(/)subhaio, and returns a distribution that is 
within 2% of the stellar profile. Thus, a rule as simple as 
selecting particles according to their potential at accretion 
is able to reproduce the radial distribution resulting from 
the complicated physics of star formation. 

Note that increasing the threshold, selects particles 
from deeper in the potential well (e.g. (j>th = 170subhaio, the 
triple dot dashed line in Fig. |4]) of the halo they were ac- 
creted in and results in a more centrally concentrated 2 = 
distribution. This is because particles with higher potentials 
at accretion are harder to strip - this requires their subhaloes 
to be on highly radial orbits such that pericentric passages 
bring them into regions where tidal forces are strong enough 
to rip them from their hosts. This occurs only towards the 
centre of the halo and results in their deposition closer to 
the halo centre. 



We have now obtained a "rule" by which DM particles 
can be selected in order to reproduce the 2 = radial dis- 
tribution of star particles. The rule is fairly simple: for each 
DM particle that is bound to a subhalo at accretion and later 
stripped, it must be sitting, at the moment of accretion, at a 
position in the potential well of its host that is deeper then 
~ 12 (jisubhaio in a radiative gas-dynamical simulation and at 
a potential that is deeper than ~ 16 ^subhaio in a DM only 
simulation. In this way a DM only simulation can be used to 
study the stellar halo without a complicated semi-analytical 
model to treat the baryons. This rule is consistent across 
our three haloes. 

Note that despite this universality, our three halos have 
very different histories. The halo of M33 had a relatively 
quiet mass accretion history, growing by smooth accretion 
and minor mergers for the past ~ 10 Gyrs. In contrast, the 
haloes of M31 and the MW experienced a more violent past 
with major mergers occurring more frequently and more re- 
cently. Our small sample size of just three haloes thus rep- 
resents a very wide variety of mass accretion histories. Al- 
though it is difficult to derive a universal relation based on 
a sample size of just three haloes, the fact that we find the 
exact same threshold value of 16(^subhaio across haloes with 
very different mass accretion history hints at the possibility 
that this is indeed a universal relation. 



3.3 The assembly of the stripped halo 

In Section 13.11 we showed that star particles and DM par- 
ticles occupy different parts of a (sub)halo's potential at 
accretion. Since star particles sit deeper in their host's po- 
tential we can infer that they are thus harder to strip and 
will thus become unbound later than DM. For each particle 
in the stripped halo, we can thus measure how many Gyrs 
after accretion the particle becomes unbound from its sub- 
structure and begins to orbit in the diffuse halo as debris. In 
Fig. [5] we show the cumulative distribution of this quantity. 
Star particles take a considerably longer amount of time to 
become unbound than their DM counterparts. For example, 
if one examines the fate of particles ~ 1 Gyr after accretion, 
we see that ~ 70% of DM particles have been stripped while 
just ~ 40% of star particles have been stripped. Some star 
particles are so deep in the potential wells of their hosts, 
that it can take up to half a Hubble time to strip them, 
while a very small fraction of DM particles that end up get- 
ting stripped will still be bound after such a long period - 
they will be stripped much earlier. 

The fact that it takes star particles significantly longer 
to be stripped than DM particles leads to the conclusion 
that the stellar halo was assembled later than the diffuse 
DM halo. In Fig. [6] we show the fraction of the stripped 
stellar and DM halo in place as function of time since the 
Big Bang. Note that after ~ 5 Gyrs, the DM halo grows 
faster that the stellar halo - reffective of the ease with which 
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Figure 5. The fraction of particles that are stripped from their subhaloes after a given period of time. In red we show star particles and 
in black, DM. The fact that at accretion star particles are deeper in their host potential is reflected by the fact that it takes longer to 
strip a given fraction of star particles than DM ones. 



DM particles are stripped. As a result, the stripped stellar 
halo is considerably younger than the DM one. In Fig. [5] 
we also show the age at which 50% of the halo debris has 
already been deposited and note that this is always earlier 
for the DM halo by around ~ 0.5 Gyrs. 

For each particle in our simulation we know the red- 
shift at which it was accreted and the redshift at which it 
was stripped. We may thus ask the question: "after accre- 
tion, what fraction of the future particle's life is spent in 
its subhalo and what fraction is spent orbiting the main 
halo as stripped debris". This question gives us a feel for 
how efficient tidal stripped is in our subhaloes. DM particles 
spend on average 10% of their post-accretion life still bound 
to their subhalo of origin and massive 90% of their post- 
accretion life orbiting as debris. As expected from the above 
arguments, stellar particles spend more time still bound to 
their subhaloes of origin: on average they spend 15% of their 
post-accretion life in subhaloes and just 85% of their time 
orbiting as debris. Note that these numbers are averages over 
all particles. Many particles - e.g. those accreted in substruc- 
tures on circular orbits - will spend more time bound to their 
substructures then orbiting as debris material. 



4 DISCUSSION AND CONCLUSIONS 

In this paper we have studied three Galaxy sized haloes 
formed in a constrained simulation of the Local Group. The 
three galaxies have approximately the same size and relative 
positions (as well as some other properties) as the observed 
Local Group members and are therefore refereed to as MW, 
M31 and M33. By construction, the galaxies are formed in 
an environment whose bulk properties (e.g. distance to a 
Virgo mass cluster) closely match observations. We use our 



constrained Local Group to focus on the similarities and 
differences in the origin of the DM and stellar halo. 

We have focused our analysis on attempting to under- 
stand why stripped DM particles orbiting in the parent halo, 
have a completely different radial distribution than ex-situ 
star particles, who were unbound by the exact same pro- 
cess. Specifically, the stars are more centrally concentrated 
then the DM, wh ich roughly follows an NFW profile (e.g. 
IWang et aLluOlOT ). Their different z = Q distribution im- 
plies that the process of tidal stripping must be operating 
on different subsets of the particles' distribution. 

We searched for a subset of the DM particles that is in- 
distinguishable from the stellar halo, in order to understand 
how these very similar components evolved differently from 
each other. We found that DM particles (in a gas-dynamical 
simulation) that had potentials at accretion of at least ~ 12 
times the potential at the subhalo's edge, are able to closely 
match the radial distribution of star particles. In a DM only 
simulation, the threshold above which selected DM particles 
reproduced the stellar halo distribution was ~ 16 times the 
potential at the subhaloes edge. 

Our threshold value is constant across the three haloes 
we have studied with very small scatter, despite the fact 
that our small sample size of three (roughly) Milky Way 
sized haloes have very different mass accretion histories. 

We have thus found a method through which the stellar 
halo may be modelled without the need of running a semi- 
analytical model to treat the baryons. By simply selecting 
those particles at accretion whose potential is greater than 
the threshold value quoted here, a set of particles can be 
identified that nearly perfectly matches the stellar distribu- 
tion. 

The implications of this work are therefore that the 
stripped stellar halo reflects the fate of material that sits 
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Figure 6. The fraction of the stripped stellar (red) and DM (black) halo, thats in place as a function of the age of the Universe. The 
dashed lines show the median of the distribution: the age at which 50% of the halo has been built up. The DM median is always earlier 
then the stellar median, reflective of the fact that at a given time, more of the DM halo is in place then the stellar halo. Half the mass 
of the stripped DM halo is in place 0.47 Gyrs, 0.44 Gyrs and 0.65 Gyrs before the stellar halo for the MW, M31 and M33 respectively. 



deep in a halo's potential well in an absolute sense, not in a 
relative one. It is important to note that selecting a fraction 
of DM particles according to their relative binding energy 
at accretion (e.g. the 10% most bound particles) will not 
successfully reproduce the z — stellar distribution unless 
the 10% most bound particles sit within a region of the po- 
tential which is greater than ~ 16 the potential at the host's 
edge. In fact selecting the 10% most bound DM particles at 
accretion returns a halo profile which follows the DM and is 
around 20% less centrally concentrated than the stars. This 
is because star formation occurs according to a local den- 
sity criteria, not according to a "global" property, such as 
a particles binding energy relative to the entire subhalo. A 
given halo will only form stars in its centre if (a) it is large 
eno ugh to retain its ga s and shield it from photo-ionization 
(e.g lBenson et al.ll2003h and (b) the potential in its centre is 
deep enough such that the gas density may trigger star for- 
mation. Thus the likelihood of star formation depends just 
on structural parameters (like concentration) and whether 
densities in a halo's centre are high enough. 

Additionally since the particles that sit in the deepest 
regions of the subhalo are harder to unbind, they will become 
unbound later - this implies that the stellar halo is younger 
then the DM halo: the star particles that compose the stellar 
halo were unbound from their substructures and deposited 
in the halo at lower z, than the DM. When we observe stars 
in the outer halo of the MW or M31, we must take care 
drawing conclusions regarding the DM haloes assembly. 

The diffuse DM halo has profoundly different proper- 
ties to the diffuse stellar halo. Its lack of central concentra- 
tion dominates its global profile. Although there are some 
detailed differences among our three haloes in the relative 
contribution from in-situ and stripped material to the stel- 



lar background, the bottom line is that DM particles can 
not serve as a proxy for star particles unless care is taken in 
their selection. 

Two main results have been presented here. The first is 
that the stars that constitute the diffuse stellar halos form 
at the bottom of the potential well of subhaloes, and hence 
are more bound at infall than the corresponding DM parti- 
cles that make the diffuse DM halo. This is a 'trivial' fact 
reproduced by all simulations of galaxy formation, yet there 
is no general consensus on the details underlying this fact. 
The novel and non-trivial other results is that one can find 
a simple mapping that enables the association of a subset of 
the DM particles with the stellar halo particles. The map- 
ping is based on one single parameter, namely the scaled 
value of the depth of the potential well within which the 
halo stars formed. Such a simple mapping of the stellar halo 
and the subset of DM particles, is a robust outcome of the 
hierarchical nature of galaxy formation. Yet, the actual value 
of the parameter that controls the mapping most probably 
varies with the particular implementation of numerical simu- 
lations of galaxy formation. This parameter might vary with 
the details of how star formation and feedback processes are 
modeled as well as on the particular numerical schemes ap- 
plied. This is posed here as an open question that we hope 
will be addressed by practitioners in the field. 
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